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Abstract 

In the class of models with small mixing between sterile and active neutrinos we place 
constraints on effective muon-to-sterile neutrino transition magnetic and electric dipole moments 
from the combined MiniBooNE results for the mass range of sterile neutrino 10 MeV < m s < 
500 MeV. Our results are valid for models with CP-violating interactions, for the Dirac and 
the Majorana sterile neutrino. Also we show that such dipole electromagnetic interactions can 
not be the main source of the anomalous events in the MiniBooNE experiment, because they 
fail to reproduce the anomalous event distribution of polar angle (meanwhile, for some values 
of magnetic and electric moments there would be a good agreement with the anomalous event 
distribution of reconstructed energy) . 

1 Introduction 

In last decades significant progress in experimental studies of neutrinos is achieved. Because of 
the incompleteness of Standard model of particle physics especially in this field, many questions 
about neutrino properties have been arisen. A Mini-Booster Neutrino Experiment (MiniBooNE) 
observes a so-called "anomaly": statistically significant excess of detected events at low energies in 
comparison with theoretical predictions [U El E] • Much efforts have been already spent and several 
hypotheses have been put forward to explain this phenomenon (see a list in e.g. [3]). 

One of the hypothesis consists in introducing sterile neutrino of mass in the range 40-600 
MeV, which is unstable with respect to radiative decay. There are two realizations of this idea 
in literature. First: sterile neutrinos are produced by flavor mixing with active (muon) neutrinos 
jH [S] at scattering of latter off nuclei due to neutral current weak interactions (un —> v s ), and 
subsequently radiatively decay due to a transition moment {y s — > uj). This idea has been already 
tested experimentally with negative result [B]. Thus, only a region of mass 400-600 MeV [7] remains 
allowed and waits for a special analysis of data from c- and b-factories [8]. Second: it was argued 
[5] that the dipole transition moment may also be responsible for both sterile neutrino production 
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and decay. In this case, it is a dipole transition moment between sterile and muon neutrinos, which 
does the main job. Anyway, the suggested explanations have renewed interest to dipole transition 
moments of hypothetical heavy sterile neutrinos. 

In this paper we extract new and improve old [TO] limits on values of sterile-to-muon neutrino 
transition dipole moments under condition that we can neglect flavor mixing between sterile and 
active neutrinos. To this end we analyze the results of the MiniBooNE detector which measures 
Charge-Current-like and Neutral-Current-like events from neutrino and antineutrino fluxes. These 
fluxes are mainly composed of muon neutrinos and antineutrinos with energies 200 MeV < E u < 
3 GeV. The detector cannot distinguish photon events from electron/positron events pQ, so neutrino 
transition dipole moments can be probed, because transition moments lead to a photoproduction 
in the detector. 

The paper is organized as follows. Section [2] describes what are neutrino transition dipole 
moments and introduces useful parametrization of them; then we calculate a cross section of active- 
to-sterile (anti) neutrino conversion on an atom and sterile (anti)neutrino radiative decay rates. 
Section [3] contains a probability function of sterile neutrino decay inside the MiniBooNE detector 
and the distributions of the expected photons of energy and polar angle. In particular, we show 
here that neutrinos, produced on nuclei and atomic electrons via exchange of massless particles 
(photons), as well as the photons from their decay are mostly forward, so fail to reproduce the 
event excess as function of polar angle [TJ [2] . Section [4] describes methods which we exploit to 
put constraints on dipole transition moments and summarizes our results and describes (possibly) 
relevant for the MiniBooNE anomaly extensions of the model. 



2 Model description 

Here we describe a model with the Dirac sterile neutrinos. The case of the Majorana neutrinos 
requires minor modifications, which are listed in the end of this Section. 

The most general Lorenz-invariant dipole electromagnetic interaction between neutrinos and an 
electromagnetic field is: 

£int = ^ E F i (A*ij + Ihdij) OpvT^Vi + h.c. , (1) 

i<j 

where are neutrino fields in the mass basis, N v is a number of neutrinos, fifj = fjLij (^dfj = dij\ 
with i j are magnetic (electric) transition dipole moments between massive neutrinos i and j, 
Ha (da) are diagonal magnetic (electric) transition dipole moments, T^ v = d^A u — d v A^ is an 
electromagnetic field strength tensor and a^ u = | (7^71/ — 7^7^), where 7 M are the Dirac matrices. 
We investigate a model with N v = 4, so we suppose an existence of one additional sterile neutrino^ 
with a mass m s . 

For our study it is convenient to work with gauge (flavor) states in the active neutrino sector. 
Sterile neutrino v s is a massive fermion neutral with respect to the Standard model gauge group. 
We are interested in a situation where flavor mixing between sterile and active neutrino is negligible, 
so that sterile neutrino in flavor basis v s and the heaviest neutrino in mass basis ^4 are almost the 
same: u± ~ v s . We denote the mass of sterile neutrino mass by m s . 

1 We extend the analysis to models with N v > 4 case in Section 4 
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Interaction ([I]) would induce two processes: an active-sterile neutrino conversion on atomic 
electrons and nuclei and sterile neutrino decay back to active neutrino (and photon). 

The transition dipole moments entering ([TJ are related to dipole moments in the flavor basis 
through relations: 

da/3 = T,i,j=l dijU^Upj , 

where U a i is a generalization of the Pontecorvo-Maki-Nakagawa-Sakata matrix: 

i 

It is convenient also to define them in another basis: 

rl ■ — V 4 rh J I* ^ ' 

In this paper we are mostly interested in the following term: 

tint = \^ {t$i + 75< 4 ) ° P xf P \ + hx. . (5) 
The undertaken analysis described bellow allows us to place constraints on the quantity: 



K 



Jk = Vl^l 2 + l<4l 2 • (6) 



Generally, electric and magnetic moments are complex quantities which have different phases. We 
parametrize them in the form: 

A = «£k cos ( a m) ex P (*^&) and d % = K % sin ( a m) ex P (ixjk) , (7) 

where A ^4, <^ r 4 and x^4 are real- valued parameters. Analogous notations are adopted for transition 
moments in the mass basis: 

H% = < cos (A£) exp and d% = k% sin (a£) exp (i X %) , (8) 

where A 4 ^, ^ and are real- valued parameters. 

In the MiniBooNE detector the conversion of muon (anti)neutrino to sterile (anti) neutrino would 
happen mostly on carbon atoms bounded in oil. In the quasi-elastic approximation the conversion 
cross section contains a form factor which depends on a single parameter t = —q 2 (where q is a 
4-momentum transfer carried by photon), and was given in |12j . It was used in previous studies of 
sterile neutrino transition magnetic moments [3 [TU] (in Ref. [U] a different nuclear part of the form 
factor was adopted for large t > 1CP 3 GeV 2 , however, the main contribution comes from atomic 
and coherent nuclear processes, where t is smaller, so the difference is insignificant). In the case 
of light nucleus (carbon: numbers of protons and neutrons are Z = 6 and A = 12 respectively) an 
analytical expression for the form factor is |12j : 

Z2ai f,i + , Za ' 4 ,f,2 , t < i ; (atomic part of the form factor) , 
' z 2 i t > to; (nuclear part of the form factor) , 

l 1+ 3 J 
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where a = 184.15 x (2.718)~ 1/2 Z" 1 / 3 / 



m e , a' = 1194 x (2.718)~ 1/2 Z- 2 / 3 /m e and d = 0.164 x 
A" 2 / 3 GeV 2 [12]; to = 7.39 x m 2 [10j (m e is electron mass). 

For the differential cross section of muon neutrino conversion into sterile neutrino of positive 
da +tU and negative cfcr_ „ chirality, we obtain: 



do-£ (£7, 0, 



ci cos 



M 4j 



4tt 



(l - sin (2A£) cos (<^ 4 



t 2 



v£ 4 (l±?;) 3 (l^cos#) , (10) 



where K^ r 4 , A^ r 4 , ^ r 4 , x^4 are defined in Eqs. ^ and ([7]), E 1 is an incident muon neutrino energy, 
v = \/E 2 — rrig/E and 9 is an angle between muon neutrino 3-momentum and sterile neutrino 



3-momentum in the laboratory frame. Equation (10) generalizes Eq. (4) of Ref. [TU] to the case of 
models where electric and magnetic transition moments are presented. 
For muon antineutrino the conversion cross section is: 



dot (E, i 



o 



(i cos 



47T 



1 + sin (2A^ r 4 ) cos ( 



G 2 (t) 
t 2 



vE 4 (l ±v) 3 (l =fcosi 



111) 



We see that neutrino and antineutrino cross sections complement each other in the following 



sense: at least one of them is non-zero for any choice of phases A^ 4 , <ffi 4 , x^ 4 if k^ 4 is non-zero. 
Combined with analogous property of sterile neutrino decay rate (see Eqs. (12) and (13)), it would 
follow that we can limit the value of k^ 4 . Note in passing that in Eqs. ( 10 ) and ( 11 ) we assume that 



quasi-elastic processes dominate and the energy of heavy neutrino is equal to the incident muon 
neutrino energy. 

Let us proceed with a description of heavy sterile (anti) neutrino radiative decay. The (anti)neutrino 
decays into some other (anti) neutrino of significantly smalleij^] mass (rrij <C m s ) and photon. For- 
mulas for differential decay rates of sterile (anti) neutrino of a given chirality (±) are: 



dT v ± 

6 7 <icos 6ry 



32vr 2 



,-m 



3 



(<) 2 (l T sin (2A£) cos (cf>% - X %) cos (0 7 )) 



(12) 



6 7 o?cos 9-y 



1 3 



(<) 2 (l ± sin (2X%) cos - cos (0 7 )) 



(13) 



where # 7 and 



■ 7 are polar and azimuthal angles of out-coming photon's 3-momentum in the sterile 
neutrino rest frame measured from the 3-momentum of heavy neutrino in the laboratory frame. As 



we see, generally there is an anisotropy in polar angle in Eqs. (12) and (13) 
The full width is (see e.g. [13]): 



i 



3 
i=l 



(14) 



The difference between the Majorana and the Dirac neutrinos just slightly modifies formulas 
and to get constraints for the Majorana case one should substitute everywhere by 2k m4 and 
by 2k 4 j. In this paper we present constraints for the Dirac neutrinos. In the case of the Majorana 
neutrinos one should divide the obtained constraints on for the Dirac case by a factor two. 



2 It is important that final (anti) neutrinos have significantly smaller mass in comparison with m s , otherwise the 
spectra get shifted towards small energies and all further calculations have to be modified, as we do in Section [4] 
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3 Theoretical predictions 



In this Section we derive approximate analytical formulas for the number of photon events associated 
with sterile neutrinos could be observed in the MiniBooNE detector. The MiniBooNE experiment 
has two operating modes (mode £ {neutrino, antineutrino}) . In the neutrino mode there is a 
dominant flux of muon neutrinos and a subdominant flux of muon antineutrinos, in the antineutrino 
mode there is a dominant flux of muon antineutrinos and a subdominant flux of muon neutrinos 
|14| (electron neutrinos don't play any role in our analysis and they are neglected in what follows). 
It should be emphasized that these two fluxes would be independent sources of photons. 

Everywhere we use following notations for angle coordinates: (0,<p) are for converted heavy neu- 
trinos in the laboratory frame, (9^,4>^) are for photons in the heavy neutrino rest frame, (Odet^det) 
are for photons in the laboratory frame. 

Bins for distributions of events of reconstructed quasi-elastic energy (evaluated from visible 
energy and polar angle, see Eq. (3) of Ref. pj)]) are: [U EJ |3]: 

bin = [E%f n , E^J e {[0.2 GeV,0.3 GeV], [0.3 GeV, 0.375 GeV], etc.} , (15) 

visible energy bins |21 E] : 

bin = [E min , E max ] e {[0.1 GeV, 0.2 GeV], [0.2 GeV, 0.3 GeV], etc.} , (16) 

and polar angle bins [TJ [2] : 

bin = [cos0 min ,cos6» maa; ] G {[-1.0,-0.8], [-0.8,-0.6], etc.} . (17) 

In each bin we compare the number of events expected due to dipole interactions with the 
number of the anomalous events (the excess is defined as difference between numbers of events 
observed and events predicted by the Standard model with three active massive neutrinos). For 
the energy range of detected photons under investigation the predictions depend monotonically 
on the energy of detected photons. The difference between visible and reconstructed quasi-elastic 
energies of photon is negative and negligible for photons with small polar angles 9d e t in comparison 
with detector resolution and bin sizes. So it is reasonable to neglect the difference between visible 
and reconstructed energies since almost all predicted events have very small polar angle Qdet (see 
Fig.§. 

The efficiency of the MiniBooNE detector e 7 (E) to register photon is determined as a piecewise 
function of visible energy [3j. So for energy distributions we set the efficiency inside each bin as 
corresponding constant. For polar angle distribution to simplify calculations we set efficiency to 
its minimal value e 7 = 7.3%, when we work on constraints (this gives a conservative constraints 
for photon energies below E 7 < 1.5 GeV). It is justified by a very moderate dependence of results 
(limits on k^) on the efficiency variation. 

To be observed sterile neutrinos has to decay inside a sphere of radius R = 5.0 m. The Mini- 
BooNE detector is located far away from the beam target, so heavy (anti)neutrino fluxes are ap- 
proximately homogeneous over the volume of spherical detector. For sterile neutrino of the energy 
E we define the average probability of the decay inside the inner sphere. Because heavy neutrinos 
could be produced everywhere inside the inner sphere, it is evaluated by averaging the probability 
for heavy neutrino to decay inside it: 

i r R r 2VR 2 -r 2 f f -xz\\ 

p(i)= w^l Mr l ('-^brJr' (18) 
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Figure 1: Schematic illustration of the conversion and decay of the v± in the MiniBooNE detector. 



where a dimensionless parameter x is defined as a product of inverse decay length and the radius 
of the inner sphere: 



S 2 



m. 



8vr ^E 2 - m? s 



R . 



Integral ( 18 ) is equal to: 



P(x) 



3 (2x 2 - 1 + (2x + 1) e~ 2x ) 
8^3 



(19) 



(20) 



Below, we obtain theoretical predictions for distributions of visible energy and polar angle and 
discuss approximations we use in calculations of constraints. 

From Eqs. (20), (10) and (12), (11) and (13) we give an analytical estimate for the number of 
events in both operating modes of the MiniBooNE experiment. We introduce two functions: the 
first corresponds to isotropic part of differential decay rates (12) and (13) (^-independent terms): 

G 2 {t) 



j^rnode ^particle 



« 4 ) 2 \ AT™? N G e, J dE cfcos - mjfC^m ^P(x)x 

x ((1 + vf (1 - cos 0) + (1 - vf (1 + cos 0)) (cos 9 cut {E max ) — cos O cu t(E m i n )) ; 



(21) 



the second corresponds to anisotropic part of differential decay rates ( 12 ) and ( 13 ) (# 7 -dependent 
term) : 



j^mode, particle 



G 2 (t) 



« 4 ) iM^Nce.jdEdcosOE^E 2 - ml(?^\E) ^P(x) x 
x ((1 + vf (1 - cos 9) - (1 - vf (1 + cos Ofj (cos 2 e cut {E max ) - cos 2 9 cut {E min fj , 



where p^^ neutrinc '(E) is a muon (anti)neutrino energy spectrum in a given operating mode |14j . 
e 7 = e 7 (0.5 (E m i n + E max )) is the detection efficiency, N~™£f e is a number of protons on target, 
Nc = (5.0/6. 1) 3 x 3.5 x 10 31 is a number of carbon atoms inside the inner sphere and a function 
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cos 9 'cut(E) is defined as: 
cos 0cut(E) 



v \m 3 



1. 



at 
at 
at 



E < 



l-v 
l+v 



2_ 



(23) 



E > f 



l+v 
l-v 



mode ,neutrino 



In a given bin the predicted number of events is defined as: 

^:t,^ ] = ( i - sin (2A M ))(^ 

+ (l + sin(2A,)) (A™±X 
where we introduced variables: 



^min i ^ m ax 

:,antineutrino 



Sill (2A 4 )Z5 rp , p 



mode,neutrino\ 
.] ) 



+ 



/n \ \ tomode,antineutrino\ 
+ S1H (2A 4 )i3 [F pi 



(24) 



sin (2A M ) = sin (2A^) cos - x^ 4 ) and sin (2A 4 ) 



3 

E 



sin(2A 4 i)cos(^-xS) • (25) 



Note that if cos 9 cut (E min ) < cos9 cut (E max ) then: 

\Qmode, particle i ^mode,particZe 



(26) 



and, consequently, the predicted number of events (see Eq. (24)) would be non-zero for all values 
of parameters A^ and A 4 . 

It is important that theoretically predicted number of events (24) depends monotonically on 
parameters k^* 4 and k\^. Provided the approximate relation: 



E«) 2 >(* 



tr V 



(27) 



i=i 



which is valid up to neglected flavor mixing between active and sterile neutrinos, we place con- 
straints from above on possible values of k 1 ^. Indeed, if we substitute everywhere the combination 

Ya=i ( k %) 2 by ( K ^ 4 ) ; it would only suppress the predicted number of events (24). Otherwise, if 
the relation (27) is invalid and hence flavor mixing between sterile and active neutrino is significant, 
constraints will depend not only on the mass of sterile neutrino but on its lifetime as well. 

In Fig. \2\ we compare the predicted number of events Af^ ode E j at parameter values m s 



50 MeV, f^ r 4 = 9.9 x 10 J fi B , t = 1/T = 1.5 x 1(T 8 s, sin(2A 4 ) = 0, sin(2A A1 ) = to observed 
energy distributions of excess events obtained in both operating modes. 

Then we investigate event distribution of polar angle. We introduce polar angle 9 de t and az- 
imuthal angle 4> de t between the muon neutrino beam axis and detected photon's 3-momentum. It 
is related to previously introduced polar angle # 7 for a given polar angle 9: 



cos I 



7 



cos 9 de t cos 6 + sin 9 det sin 9 cos § det - v 
1 - v (cos 9 det cos 9 + sin 9 de t sin 9 cos cf> det ) 



(28) 



To define the predicted number of events we need to cast Eqs. ( 12 ) and ( 13 ) in terms of new angles 
(9 det ,4> det ). The integration measure is transformed as: 



efcos 9 



(1 



cfcos 9 det d4> det 



(l-v (cos 9 det cos 9 + sin 9 det sin 9 cos 4> de t)Y 



(29) 
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Figure 2: Energy distributions of excess events obtained in neutrino (left) and antineutrino (right) with 
error bars |3] and corresponding predicted spectra of photons Af^ ode E i from sterile neutrino decays at 
parameter values m s = 50 MeV, «£ r 4 = 9.9 x l(T 9 /iB, r=l/T= lXx 10" 8 s, sin (2A 4 ) = 0, sin (2A M ) = 0. 

Using Eqs. (28) and (29) we transform decay rates (12) and (13) into new variables 9 de t an d 4>det- 
Analogously to the previous case we introduce two functions: 

?J mode,particle / tr \ 2 ^_ irmode *r f (l - V 2 ) dE rfcOS 9 dcOS 9 det d(j) dct 

[cos9 — ' cose — 1 " W 87/ C J (l- v (cos 9 det cos + sin 0^ sin 9 cos <^ et )) 2 X 

x&y/EP-mfatfZSf'iE) ^ ((1 + ^) 3 (1 - cos 9) + (1 - V ) 3 (1 + cos 9)) P (x) x (30) 

x@(l + v (cos det cos + sin 9 det sin cos <j) det ) - 2 ^ /l u E rrvm j ^ 
\ m s J 

ymode,particle = (tr \ 2 ^kfmodejq f (l - f 2 ) dcOS rfCQS 9 det d4> det ^ 

[c OS e min , CO B0 max ] ~ { K ^) 87r p.o.t. cj (cog — - - — — — - 

x (cos 9 det cos + sin 9 det sin cos <p det - v) E 3 ^Je 2 - m 2 e 1 p p ^ r t d % f e {E) G ^ x 

(31) 

x ((1 + v) 3 (1 - cos 0) - (1 - vf (1 + cos 0)) P (x) x 

x9 ( 1 + v (cos 9 det cos + sin 9 det sin cos <?W) - 2 ^ /l ^ j 
\ m s J 

where we integrate over <fi from to 27r and over cos 9 de t from cos 9 max to cos 9 m i n , O is the step 
function, which selects events with visible photon energies E^ > E m i n = 0.2 GeV. Then theoretical 
prediction for a number of events by cos 9 de t with visible photon energies E^ > E m i n is: 

\fmode (i „:„ (n \ \\ fj ,mode, particle _ , tr> \ \^,mode, particle \, 

* V fcos0 min ,cos0 mM; ] - V 1 sm K^n)) ^ M [co S roin ,cos0 mtra ] S111 V ZA 4j ^[cosfl^costW 1 + 

+ (1 + sin (2A M )) ( U ™^ article ft . + sin (2\AVr de ft > partlde ft 



,COS m ax 
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Figure [3] demonstrates that almost all photons would be produced in the forward direction. To 
properly constrain from polar angle distributions one needs to know the uncertainties for the 
excess events. However, we did not find these estimates in literature, and thus do not use polar 
angle distribution to place limit on k^. Note in passing that the excess distribution of cos 6 det 
remains intact with decrease of sterile neutrino lifetime (e.g. if other decay modes are introduced). 
Indeed, the decrease of lifetime would just add more energetic photons and increase the peak in 



the forward direction, because the function (20) suppresses high energetic particles stronger than 



low energetic ones and has a limit equal to one at big arguments. 
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Figure 3: Theoretically predicted events as a function of cos 9 det in the neutrino operating mode at 



50 MeV, 



1 x 10" 



i - ' <<> Mb, t = 1/r = 1.5 x 1CT 6 , sin(2A M ) = 0, sin(2A 4 ) = (solid line) and the 
observed excess (the dashed line) [3] . 



4 Analysis methods and results 



As follows from fitting theoretically predicted numbers of events (32) to the excess distributions, 
the anomalies [TJ [3] cannot be entirely explained only by the effective dipole electromagnetic 
interaction ([5j). For this reason we: (i) work on constraints, but not on the expected values of 
parameters, which will explain anomalies; (ii) analyze each bin and each mode independently. 
We find numbers of events as functions of /O, sinA^ and sinA4 for all bins and modes from 



(24). After these calculations, using known uncertainties of the experimental results for each bin 
and mode, we investigate the predicted numbers of events as functions of sinA^, sinA4 and m s . 
In each bin we constrain «^ r 4 (sin A^, sin A4, m s ) at 95% C.L.. Then we finally constrain for a 
given m s by the minimum value over all bins and modes and the maximum value over continuous 
variables sin A M and sin A4 of the constraints on (sin A M , sin A4, m s ). 

In Fig. |i| (left) the upper limit on k^ a at 95% C.L. is shown and in Fig. |i] (right) we present con- 



straint on sterile neutrino mean life time, derived from the upper limit on fO (provided Eq. (27)). 
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We see that data from the MiniBooNE detector allow to improve constraints [TO] for the mass range 
m s < 350 MeV. 




Figure 4: In the left plot the shaded area in (m s , k*^ 4 ) is excluded at 95% C.L. from our analysis. The 
region above the dashed line was excluded in Ref. [TU]. In the right plot we present constraints on sterile 
neutrino lifetime obtained from constraints on k^. 



Now let us discuss a modification of the previous analysis to the models with more sterile 
neutrinos (N u > 5). In particular v s may decay into 7 and new sterile neutrino(s) u 8 a (j > 5), if 
kinematically allowed, i.e. m s > m s j. Eqs. (12) and (13) take different forms: 



dr 



± 



1 



6 7 dcos # 7 



+ 



32vr 2 
1 



in. 



2 

J2 ( K ti) (l =F sin (2X%) cos (<f>% - xtsj) cos 1 



i=l 

N v 



~m Q 

32^ 2 s 



mi 



m„ — 



3=5 



mi 



x (l =F sin (2A| r Sii ) cos (^>% ;j - xt,j) cos (ft 
' l],± ' " 3 E(4) 2 (l±sin(2AS)cos(^-X^) cosl 



defied cos 9~, 32tt 2 



1 



nr 



nni" 

32vr 2 s 



i=l 
i=5 



m, — m 



m 



2 ) ( K 4s,j) 



1 ± sin (2\ts,j) cos (4>t s j ~ xts,j) cos ' 



+ 



) + 



The energy of the outgoing photon as the function of E and cos 9 is determined through: 



E-, 



nr 



s,3 



mi 



1 + v cos 9 



(33) 



(34) 



(35) 
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And consequently we have to introduce modifications of the cutoff function 9 cu t (see Eq. (23)) 
which now depends on m s j: 



( 



cos 9^ t (E) = cos 9, 



cut 



E 



m . 
1 s,J 
J- 9 



(36) 



We are ready now to discuss the second mentioned in Introduction scenario [9] suggested to 
explain MiniBooNE and LSND anomalies. Initially it contained flavor mixing between sterile and 
muon neutrinos, but such mixing was found inconsistent with kaon decays [6]. Then the authors 
of [5] have introduced two sterile Majorana neutrinos: m s = 50 MeV, 5 MeV < m s ^ < 10 MeV and 
two transition magnetic dipole moments /x^ r 4 = 2.4 x lO -9 ^ and ^45 = 2.4 x 10~ 8 /i#, the latter 
saturates the heavy sterile neutrino decay rate via v s — > v s jl- Equation (35) implies that massless 



approximation is reasonable for such hierarchy of m s and m s ^. To test the case with two sterile 
neutrinos we fix the proposed relation between magnetic moments //* r 4 = 0.1 x /X45. And using 



Eq. (24) we place constraint /4s < 2.1 x 10 _8 /xb (95% C.L.) for the Majorana neutrinos, which 
is marginally consistent with the given value ^45. However, evidently such model is disfavored, 
because additional source of photons would give a peak in the forward direction. Extensions of the 
obtained constrains to the models with several neutrinos are straightforward. More careful analysis 
is required to test models with almost degenerate (in masses) sterile neutrinos. In particular, in 



K^^j most 



the model with two degenerate neutrinos ~ m s ^ for (k%) 2 (l — rn 2 s ^/m 2 ^j 3> ( 
photons from decays —> v s ,bl could have energies comparable or below lower energy cut off 
at 200 MeV. In such specific case sterile neutrino could avoid ours constraints for k% and k^. 
But such situation is excluded as the explanation of the excess, because experimental data shows 
smooth angular dependence |3 . 



5 Conclusion 

Combined analysis from different focusing regimes of the MiniBooNE detector have allowed us to 
investigate models with heavy sterile neutrinos. Particularly we renewed and generalized constraints 
on electromagnetic transition dipole moments. The main results are presented in Fig.|4j The 
obtained analytical formulas are rather general and applicable for other neutrinos experiments. 

We have shown that sterile neutrino of mass 10 MeV < m s < 500 MeV with transition dipole 
moments and without significant mixing with active neutrinos of the mass in the region 10 MeV < 
m s < 500 MeV, cannot explain anomalies if one regards MiniBooNE polar angle distribution of the 
anomalous events. 
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